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Heavy quark bound state systems, mesons and hybrids, are discussed in the frame- 
work of the QCD effective field theory called potential Non-Relativistic QCD. 



1 The physical systems 

Among all the hadrons, the ones that should be the simplest to analyze are 
those entirely composed by heavy quarks (and gluons). As it is evident from 
the spectra (e.g. by comparing the characteristic energy levels splitting with 
the value of the mass of the quarks in heavy quarkonium), these systems are 
nonrelativistic. Thus, they can be described in first approximation using a 
Schrodinger equation with a potential interaction. This amounts to saying 
that the heavy quark bound state is characterized by three energy scales, hi- 
erarchically ordered by the quark velocity v ^ 1: the mass m (hard scale), 
the relative momentum mv (soft scale) and the binding energy mv^ (ultrasoft 
scale). These scales typically get mixed in any Feynman diagram bound state 
calculation and originate technical complications (the same happens in QED 
e.g. for positronium). On the other hand, even in a lattice calculation, it is 
difficult to handle physical systems that possess two very different scales, say 
Q ^ q, since a really demanding relation should then hold for the lattice size 
L and the lattice step a: <^ q <^ Q <^ . 

In QCD, a further conceptual complication arises due to the existence of 
a nonperturbative scale, call it Aqcd, the scale at which the nonperturbative 
effects become dominant. For heavy quarks only the hard scale m is surely 
bigger than Aqcd and can be treated perturbatively. Many attempts have been 
made to properly consider the effect of the nonperturbative dynamics on the 
heavy quark energy levels. On the one hand, phenomenological or QCD derived 
confining potential models have been used inside a Schrodinger equation. Such 
a picture suffers from many ambiguities, it is strongly model dependent, it 
consists in a by hand superposition of perturbative and nonperturbative effects, 
and it leaves out the physics connected to the retardation effects, that's to say 
the physics at the ultrasoft (US) scale. On the other hand, the contribution to 
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e Coulombic energy levels due to local condensates has also been evaluated 
This correction is of nonpotential type (it is analogue to the Lamb shift 
effect) and it grows out of control for the excited levels. 

Summarizing, the existence of different scales and the nontrivial features 
related to perturbative (hard scale) and nonperturbative (low scale) effects 
and to potential and nonpotential contributions, complicate considerably the 
description of the heavy quark bound states. I will show here how it is pos- 
sible to take advantage of the existence of this hierarchy of widely separated 
energy scales to construct QCD effective field theories (EFT) with less and less 
degrees of freedom. This leads ultimately to a field theory derived quantum 
mechanical description of these systems. The corresponding EFT is called pN- 
RQCD (potential NonRelativistic QCD)a'El. Here, all the dynamical regimes 
are organized in a systematic expansion. 

2 The QCD effective theories for these systems 

The idea is that since the typical scales of a nonrelativistic bound state system 
are widely separated, it is possible to perform an expansion of one scale in 
terms of the others. Roughly speaking, first one expands in 1/m and then one 
performs aJx.ppapansion in the inverse of the soft scale, the so called multipole 
expansion otla'El. The effective theory supplies us with the procedure to make 
this expansion consistent with the ultraviolet behaviour of QCD and consistent 
with a systematic power counting in the small expansion parameter v. As I 
will explain below, it depends on the physical system in consideration and on 
the actual position of Aqcd into this hierarchy of scales, if it is possible or not 
to take advantage of the further simplification of performing the calculations 
in a perturbative expansion in as- 

2.1 NRQCD 

First we pass from QCD to NRQCD by integrating out the hard scale m. This 
involves performing an expansion in 1/m, which, in the two-fermions sector is 
of the type of the Foldy-Wouthuysen transformation. Since we are modifying 
the ultraviolet behaviour of the theory, matching coefficients and new operators 
have to be added in order to mock up the effects of heavy particles and high 
energy modes into the low energy EFT. In this case, since the scale of the mass 
of the heavy quark is perturbative, the scale /i of the matching from QCD to 
NRQCD, mv < ^ < m, lies also in the perturbative regime. Then, the hard 
scale is integrated out by comparing on shell amplitudes, order by order in 1/m 
and in as, in QCD and in NRQCD. The difference is encoded into the matching 
coefficients that typically depend non-analytically on the scale which has been 
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integrated out, in this case to. We work here in dimensional regularization, 
MS schema, and with quark pole masses. Up to order 1/to^ the Lagrangian 
of NRQCDB reads: 

^ ./^ D2 D4 SB DEED 

£ = [iDa + C2- h C4— ^ + cpg h cog 5-^ (1) 

\ 2to oto'^ to om^ 

S • (D X E - E X D)\ , . , 

+icsg ip + antiquark terms + terms with light quarks 

4to^ / 

— !!lp<^ pa^'^ I n2 pa/iiy , /■ pa pb pc 

^ ^ fj,u^ ^ ^2 A"' ^2"''"''^ A"' ''a 

TO^ TO^ TO^ TO^ 

■0 and X being respectively the quark and the antiquark field; D^^ is the covari- 
ant derivative, E and B are chromoelectric and chromomagnetic fields, F'^'^ 
is the gluon field strength, S is the total spin, T° is the color generator, g is 
the coupling constant. Ci,bi,di are matching coefficients. They depend on /i 
and TO and are known in the literature at a different level of precision. The fi 
dependence in the matching coefficients cancels against the fi dependence of 
the operators in the Lagrangian. The gluonic part in the third line comes from 
the vacuum polarization while the last line contains four quark operators. 

Many applications of NRQGD have been made both in lattice calculations 
and in the continuum, see e.g. Ell3. Here, however I just introduced NRQCD 
as a step towards pNRQCD. I recall in fact that in NRQCD two dynamical 
scales are still present, mv and tow^. These scales get entangled and obscure 
the power counting, i.e. the matrix elements of the operators in (|^) do not 
have a unique power counting but they also contribute to subleading orders 
in V. In the next section I will show how it is possible to integrate out also 
the momentum scale, obtaining an effective theory at the ultrasoft scale. This 
effective theory is called pNRQCD. 



2.2 pNRQCD 

When we integrate out the soft scale mv, there are two possible situations 
depending on the relative position of Aqcd with respect to the other scales. If 
Aqcd is smaller than mv, then also the matching from NRQCD to pNRQCD 
may be performed in perturbation theory. If Aqcd — "tiv the soft scale is 
already nonperturbative and the matching to pNRQCD has to be performed 
avoiding any expansion in a^. Roughly speaking, we can say that the lowest 
excitations of quarkonium belong to the first situation while the excited states 
belong to the second situation. Indeed, the typical radius of the bound state 
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is proportional to the inverse of the soft scale r ~ l/mv and thus for the 
lowest states the condition mv ^ Aqcd may be fulfilled. In Fig.l I present 
the typical radius of various mesons against the curve of a phenomenological 
potential that displays a nonperturbative (linear) behaviour around 0.2 fm. In 
the following I will discuss both situations. 

3 ■ 
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Figure 1: rqq of various mesons vs the Cornell potential, from^. 



3 pNRQCD for mv > Aqcd 

We denote by R = (xi + X2)/2 the center of mass of the QQ system and by 
r = Xi — X2 the relative distance. The matching to pNRQCD is perturbative. 
At the scale of the matching ^' {mv /x' » mw^, Aqcd) we have still quarks 
and gluons. The effective degrees of freedom are: QQ states (that can be 
decomposed into a singlet S'(R, r, t) and an octet 0(R, r, t) under color trans- 
formations) with energy of order of the next relevant scale, 0{Aqcd-, mv"^) and 
momentum^ p of order 0{mv), plus ultrasoft gluons A^(R, t) with energy and 
momentum of order 0(Aqcd, m-^^^)- Notice that all the gluon fields are mul- 
tipole expanded. The Lagrangian is then an expansion in the small quantities 
p/m, 1/rm and in 0(Aqcd, wi'^) x r. 

3.1 The pNRQCD Lagrangian 

At thepjj^xt-to-leading order (NLO) in the multipole expansion, i.e. at 0(r), 
we gett 
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4nrqcd = -iKuF''"' + Tr <! St ( - ^ - - y: ^ ) S ^ (2) 



m ' — ' m" 

n—l 



''Notice that although, for simplicity, we described the matching between NRQCD and pN- 
RQCD as integrating out the soft scale, it should be clear that the relative momentum p of 
the quarks is still soft. 



+gVATr {Oh • E S + S^r • E O} + ^-yTr {OV ■ E O + O^Or ■ E} . 

The matching coefficients Vj are functions of /i, fi' , r, p, Si, S2. The equivalence 
of pNRQCD to NRQCD, and hence to QCD, is enforced by requiring the 
Green functions of both effective theories to be equal (matching). In practice, 
appropriate off shell amplitudes are compared in NRQCD and in pNRQCD, 
order by order in the expansion in 1/m, as and in the multipole expansion. The 
difference is encoded in these potential-like matching coefficients that depend 
non-analytically on the scale that has been integrated out (in this case r). 
Recalling that r ~ 1/mv and that the operators count like the next relevant 
scale, 0{mv^ , Aqcb), to the power of the dimension, it follows that each term 
in the pNRQCD Lagrangian has a definite power counting. This feature makes 
LpNRQCD the most suitable tool for a bound state calculation: being interested 
in knowing the energy levels up to some power w", we just need to evaluate 
the contributions of this size in the Lagrangian. The singlet sector of the 
Lagrangian would give rise to equations of motion of the Schrodinger type, 
while the last line of (||) contains retardation (or nonpotential) effects that start 
at the NLO in the multipole expansion. At this order the nonpotential effects 
come from the singlet-octet interaction mediated by a ultrasoft chromoelectric 
field. 

Summarizing, pNRQCD is equivalent to QCD, it has potential terms, 
thus embracing potential models, and it has ultrasoft gluons incorporated in 
a second-quantized, gauge- invariant and systematic way. Moreover, the power 
counting in v is explicit and subsequent corrections both in the 1/m expan- 
sion as in the multipole expansion can be added systematically. Hard and soft 
effects are separated from the ultrasoft (low-energy) effects. 

The Feynman rules corresponding to the Lagrangian (|^) are shown in Fig. 2. 

At zero order in the multipole expansion the singlet and the octet decou- 
ple and their equations of motion turn out to be Schrodinger-like. Then, the 
first question comes: is the singlet matching coefficient Vs equal to the static 
heavy quark potential? The answer depends on the actual ratio of mv'^ and 
Aqcd- In the EFT language the potential is defined upon the integration 
of all the scales up to the ultrasoft scale mv^. We can imagine two differ- 
ent situations: if Aqcd ^ mv^, then Vs is the static heavy quark potential; 
if mv <C Aqcd ^ mv^, then we have to integrate out also Aqcd in order 
to get the potential. In this integration the potential acquires short range 
nonperturbative contributions as I will show in Sec. 3. 3. 
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A second interesting question is what is the relation between Vg and the 
energy Eg between quark static sources, which is defined as 



^.(r-)= lim -ln{Wu), 

1 — *CO ± 



(3) 



being Wu the static Wilson loop of size r x T, and the symbol ( ) being the 
average over the gauge fields. Eg is often used as a static potential inside the 
Schrodingcr equation, assuming that the Born-Oppcnhcimcr approximation 
holds. We answer these questions by performing explicitly the singlet matching 
at order l/wP and at the NLO in the multipole expansion. 



I 
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singlet propagator 
octet propagator 
singlet-octet vertex 
octet-octet vertex 
Coulomb octet-octet vertex 



Figure 2: Propagators and vertices of the pNRQCD Lagrangian (at order l/m*^ and at the 
NLO in the multipole expansion). 



3.2 The matching procedure: the singlet fAqcD < mv'^) 

The matching can be done once the interpolating fields for S and O" have been 
identified in NRQCD. The former need to have the same quantum numbers and 
the same transformation properties as the latter. The correspondence is not 
one-to-one. Given an interpolating field in NRQCD, there is an infinite number 
of combinations of singlet and octet wave- functions with ultrasoft fields, which 
have the same quantum numbers and, therefore, have a non vanishing overlap 
with the NRQCD operator. However, the operators in pNRQCD can be orga- 
nized according to the counting of the multipole expansion. For instance, for 
the singlet we have 

xHx2, i)<^(x2, xi; t)V(xi, = Zll\r)S{n, r, t)+Z]ll(r)r r.E«(R, t)0\B., r, t)+. . . , 

(4) 
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and for the octet 

X^(x2,i)0(x2,R;OT>(R,xi;i)V'(xi,<) = Zy'{r)0'^{R,r,t) (5) 

+ z'J'^{r)rr-E-{R,t)S{R,r,t) + ..., 

(j)(y, x; i) = P exp{ig J^ds (y — x) • A(x — s(x — y), t)}, Zi being normalization 
factors. These operators guarantee a leading overlap with the singlet and the 
octet wave-functions respectively. Higher order corrections are suppressed in 
the multipole expansion. The expressions for the octet can be made manifestly 
gauge-invariant by inserting a magnetic field in place of the pure color matrix. 
The fact that the matching can be done in a completely gauge-invariant way 
enables us to generalize pNRQCD to the case in which Aqcd — rnv, i.e. to 
the nonperturbative matching, cf. Sec. 5. 

Now, in order to get the singlet potential, we compare 4-quark Green 
functions. From (^, we take in NRQCD 

/ = (53(xi-yi)j3(x2-y2)(W^n), (6) 

and we equate to the singlet propagator in pNRQCD at NLO in the mul- 
tipole expansion (cf. Fig. 3 for a diagrammatic representation) 

/ = Z,(r)<53(xi - yi)<53(x2 - y2)e-^^^=W x (7) 

^ J^c J-T/2 J-T/2 J 

where (j)'^'^^ is a Schwinger (straight-line) string in the adjoint representation 
and fields with only temporal argument are evaluated in the centre-of-mass 
coordinate. From here one gets the singlet static potential Vs (and the singlet 
wave- function normalization Zg)'. 

Vs (r) = (r) I 2_ioop+NNNLL + W T2" ' 

where Es has been defined in Eq. (^)^ We note that Vg and Eg would coincide 
in QED and that therefore the effect we are studying here is a genuine QCD 
feature. The leading log contributions to the NNNLO arise from three loops 

'^The non-relativistic limit described by the Schrodinger equation with the static potential 
is called 'leading-order' (LO); contributions corresponding to corrections of order v" to this 
limit are called N"LO. LL means 'leading-log'. In the perturbative situation v = as- 
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NRQCD pNRQCD 

Figure 3: The matching between the Wilson loop in NRQCD and the singlet propagator in 
pNRQCD. 

diagrams that are infrared divergent. The final result reads 

Vs{r,^j.) = -Cf , 9 

r 

av{r, ^J.') = as{r) 1 + aias(r) + a2{as{r))'^ + — In/xV 

L TT iZ 

This is the LL three loop evaluation of the static heavy quark potential in the 
case Aqcd < mv"^. The one loop and two-loop contributions come fromB and 
the three loop LL fromi. 

This situation is expected to hold for toponium and for the bottomonium 
(charmonium?) ground state. The explicit n' dependence of Vs originates from 
the fact that the US degrees of freedom (which have the same scale of the ki- 
netic energy and therefore do not belong to the potential) have been explicitly 
subtracted out from the static Wilson loop. This fact is not surprising if we 
understand the heavy quarkonium potential as a matching coefficient of pN- 
RQCD. As a consequence even in a purely perturbative regime the static heavy 
quarkonium potential (as well as ay) turns out to be an infrared sensitive quan- 
tity. In this situation nonperturbative effects are only of non-potential nature 
and appear in the form of local (a la Leutwyler-Voloshin) or nonlocal conden- 
sates. When calculating any physical observable, the fi' dependence canpds 
against //'-dependent contributions coming from the dynamical US gluonsEJ. 

In a similar way, by matching appropriate Green functions, [We, can get the 
octet matching potential Vo and the other matching potentials cllly. 
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3.3 The singlet potential in the situation mv ^ Aqcd ^ 

Since in this situation there is a physical scale (Aqcd) above the US scale, a 
potential can be properly defined only once this scale has been integrated out. 
At the NLO in the multipole expansion we get 

V,ir) = -CF^^^^-^^TFV2ir)J|lite'^'^^^^^ 

(10) 

Therefore, the heavy quarkonium static potential Vg is given in this situation by 
the sum of the purely perturbative piece calculated in Eq. (||) and a new term 



8 



carrying also nonperturbative contributions (contained into non-local gluon 
field correlators). This last one can be organized as a series of power of r" 
by expanding exp{~it{Vo - Vs)} = 1 - it{Vo - Vs) + . . . (since t ~ 1/Aqcd, 
Vq — Vs ~ mv^). Typically the nonperturbative piece of Eq. ( |l0| ) absorbs the ^' 
dependence of ay so that the resulting potential Vs is now scale independent. 
We notice that the leading nonperturbative term could be as important as 
the perturbative potential once the power counting is established and, if so, it 
should be kept exact when solving the Schrodinger equation. In Table 1 we 
summarize the different kinematic situations. 



mv 




potential 


ultrasoft corrections 


» Aqcd 
» Aqcd 
> Aqcd 


> Aqcd 
^ Aqcd 
< Aqcd 


perturbative 
perturbative 
perturbative + 
short-range nonpert. 


local condensates 
non-local condensates 
No US (if light quarks 
are not considered) 



Tabic 1: Summary of the different kinematic situations. 



We can also consider Eq.([l0|) in the limit of static sources. In this case 
there are no dynamical scales involved and therefore in the short range we have 
exp{-it{Vo~Vs)]{^%t)(t>{t,Qf^''E\0)) ~ exp{^it{Vo~Vs)}{E{Of). This 
gives 

Es{r) ^ -C,^^^^ - '1^{E\Q)) (11) 
r 3 N^a^ 



which coincides with 



4 Application of pNRQCD (Aqcd < rnv) 

• Quarkonium spectrum at leading-log of NNNLO 

The complete leadingJcg terms of the NNNLO corrections to pNRQCD 
have been calculated EEl. As a byproduct the leading logs at 0(r7ia^) in 
the heavy quarkonium spectrum have been obtained. When Aqcd ^ 
TOQig, these leading logs give the complete 0(ma^ In a^) corrections to 
the heavy quarkonium spectrunL. (plus the nonpotential contributions 
coming from local condensates) 113. The result is important at least for 
ti production and T physics. In the first case this result is a first step 
towards the goal of reaching a 100 MeV sensitivity on the top quark 
mass from the ti cross section near threshold to be measured at future 
linear colliders. In the second case it will improve our knowledge on the 
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b mass. In botlL.|Ca,ses the LL contributions of NNNLO have been found 
to be relevant BE3. 

• Renormalons 

The infrared sensitivity of the static potential can also be expressed in 
the renormalons language, i.e. we can say that Vs, as defined in Eq.(^, 
suffers from IR renormalons ambiguities with the following structure 

Vs{r)\ IR rcn — Co + Car^ + . . . (12) 

The constant Co ~ Aqcd is known to be cancelled by the pole mass 
IR renormalon (2mpoio|iRrcn — ~Co). While Eq. (j^j provides us with 
the explicit expression for the operator which absorbs the C2 ^ ^qcd 
ambiguityi. The renormalon cancellation issue plays an important role 
in quarkonium phenomenology cJ. 

• Hybrids and gluelumps 

pNRQCD gives model independent predictions on the behaviour of the 
hybrid static potentials. In particular it predicts for these potentials 
an octet behaviour at very short distances and—it correctly states all 
the degeneration patterns in the small r region 0. Moreover, it allows 
to relate the mass of the gluelumps to th£_.cprrelation lengths of some 
nonlocal vacuum field strength correlatorsEJa. 

• Quarkonium scattering, Van Der Waals forces, Quarkonium Production. 
Work is in progress on these applications. 

• Renormalization group improvement _ 
Up to now, it has been performed at the level of the static potential E3. 



5 pNRQCD for A ~ mw 

In this case the potential interaction is dominated by nonperturbative effects. 
This is the most interesting situation, in which most of the mesons lie (cf. 
Fig.l). 

A large effort has been made in the last decades in order to ahtaiii from 
QCD the nonperturbative potentials in the Wilson loop approachtOHB. pN- 
RQCD allows us to. obtain via a nonperturbative matching all the nonpertur- 
bative potentials ta. From this respect I want to discuss few concepts and 
results. 

In pNRQCD a potential picture for heavy quarkonium emerges at the lead- 
ing order in the US expansion under the condition that the matching between 
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NRQCD and pNRQCD can be performed within an expansion in 1/m. The 
gluonic excitations (hybrids and gluebaUs) that form a gap of order Aqcd with 
respect to the quarkonium state can be integrated out and the potentials foUow 
in an unambiguous and systematic way from the nonperturbative matching to 
pNRQCD. Thus, we recover the quark model from pNRQCD. The US degrees 
of freedom in this case are not coloured gluons but US gluonic excitations be- 
tween heavy quarks and pions. They can be systematically included and may 
eventually affect the leading potential picture. Let us consider for instance 
the singlet matching potential. Disregarding the US corrections, we have the 
identification 

Vs{r) ^ Yyu,^^\n{W(r XT)). 

US corrections to this formula are due to pions and US gluonic excitations 
between heavy quarks, and may be included in the same way as the effects 
due to US gluons have been included in the perturbative situation through Eq. 

®- n 

The complete 1 /m? potential has been calculated along these lines 113. 
There are many appealing and interesting features of this procedure. The 
matching calculation is performed in the way of the quantum mechanical per- 
turbation theory on the QCD Hamiltonian (where perturbations are counted 
in orders of 1/to and not of Us) and only at a later stage the relation with the 
Wilson loop and field insertions is established. This allows us to have a control 
on the Fock states of the problem and on the contributions coming from glu- 
onic excitations in the intermediate states. At the end, all the expressions are 
again given in terms of Wilson loops, which can be evaluated on the lattice or 
in QCD vacuum models E3. The potentials turn out to be naturally factorized 
in a hard part (the matching coefficients at the hard scales inherited by pN- 
RQCD from NRQCD) and a low energy part (the Wilson loops expressions). 
The power counting may turn out to be quite different from the perturbative 
(QED-like) situation. 

6 Conclusion and outlook 

I have shown that it is possible to construct systematically and within a con- 
trolled expansion an effective theory of QCD, which describes heavy quark 
bound states. All known perturbative and nonperturbative regimes (poten- 
tial, nonpotential), are dynamically present in the theory, which is equivalent 
to QCD. I have presented many applications of pNRQCD in the situation 
Aqcd Tiw. In the situation Aqcd — t^v, I have discussed how pNRQCD 
allows us to systematically factorize the nonperturbative heavy quark dynam- 
ics. 
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